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Optimization of Excitation Waveform for Nonlinear 
Transmit-Receive Systems 

Cross-Reference to Related Applications 

[0001] This application is based upon and claims priority to U.S. Provisional 
Application Serial No. 60/397,378, entitled "Optimization of the Excitation Waveform 
for Transmit-Receive Systems," filed July 19, 2002, the entire content of which is 
incorporated herein by reference. 

Background of Invention 
[0002] Field of Invention 

[0003] This invention relates to transmit-receive systems and, more particularly, to 
optimizing the excitation waveform that is used in such systems to maximize the 
output for a given constraint on the input. 

[0004] Description of Related Art 

[0005] Systems employing a transmitter and receiver are widely used. Applications 
include telemetry, imaging and communications. A specific example of an imaging 
application is ultrasonic transmission-mode tomography used in breast 
mammography. 

[0006] A signal of a selected waveform is typically delivered to the transmitter, 
emitted by the transmitter, propagated though a medium, and received at the receiver. 
Information of interest is then extracted. 

[0007] To maximize the signal-to-noise ratio, it is often desirable to maximize the 
output at the receiver based on a given constraint on the input to the transmitter. 

[0008] One approach has been to optimize the design of the individual hardware 
components, such as the transducers and electronics. For example, the thickness and 
impedance of the various matching layers in the transducers have been adjusted. 
Appropriate tuning circuitry has also been introduced to maximize the effect of the 
excitation waveform on the transmitted transducer. 



LAS99 1299341-1.064693.0074 



-1- 



► 1 



Matter No. 64693-074 



[0009] These types of tuning methods can be problematic. The impedance of the 
electronic circuits associated with the transmitter and receiver may be different, 
requiring different tuning conditions during transmission and reception. Precise 
knowledge of the dynamic characteristics of the transducers may also be needed in 
conjunction with the associated electronics. 

[0010] Rather than physically tuning system components, transmit-receive systems 
have been modeled as linear systems. Based on this modeling, the center frequency, 
shape, duration and amplitude of the excitation waveform has been adjusted to 
maximize the peak received value predicted by these models. 

[0011] Unfortunately, there are often inherent nonlinearities in transmit-receive 
systems, including in ultrasonic transmit-receive systems, that can cause the results of 
this linear modeling approach to be less than optimal. In turn, the greatest possible 
signal-to-noise ratio may not be achieved. 

Summary of Invention 

[0012] A process for optimizing the excitation waveform that is delivered to an 
ultrasonic transmitter that, together with an ultrasonic receiver, form part of a nonlinear 
ultrasonic transmission and reception system. A transmission test signal may be 
delivered to the ultrasonic transmitter. A received test signal may be generated from 
the ultrasonic receiver that is a nonlinear function of the transmission test signal. A 
nonlinear model of the nonlinear function may be developed from the transmission test 
signal and the received test signal. An optimum excitation signal for the ultrasonic 
transmitter may be determined that substantially maximizes the signal generated by 
the ultrasonic receiver based on the model and based on a specified constraint on the 
excitation signal. 

[0013] The kernel functions of the nonlinear function may be determined as part of 
the process of developing the model. An algorithm may be used in determining the 
kernel functions. Principal dynamic modes of the nonlinear function may also be 
determined as part of this process and may be based on the kernel functions. 
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[0014] The time inversion of one or more of the principal dynamic modes may be 
determined as part of determining the optimum excitation signal. 

[0015] One or more of the kernel functions may be excluded when calculating the 
time inversion. 

[0016] A Laguerre-Volterra network may be used in developing the nonlinear model. 
Parameters of the Laguerre-Volterra network may be adjusted to minimize the mean- 
squared error between the signal predicted by the network and the received test 
signal. The adjustment may be an iterative process. 

[0017] The nonlinear model may include a linear filter followed by a static 
nonlinearity. 

[0018] The specified constraint on the excitation signal may include a constraint on 
the amplitude or power of the excitation signal. 

[0019] The amplitude or power of the signal generated by the ultrasonic receiver may 
be maximized as part of determining the optimum excitation signal. 

[0020] The transmission test signal may be a wideband signal. The wideband signal 
may cover the bandwidth over which the ultrasonic transmitter is configured to 
operate. The wideband signal may be white noise or a chirp. 

[0021] The transmission test signal may cover the dynamic range over which the 
ultrasonic transmitter is configured to operate. 

[0022] The ultrasonic transmitter may be excited with the optimum excitation signal. 
A breast may be placed between the ultrasonic transmitter and the ultrasonic receiver 
during this process. The signal received by the ultrasonic receiver in response to the 
optimum excitation signal may be analyzed to create an image of the breast. 

[0023] An ultrasonic imaging system for generating an image of tissue may include 
an ultrasonic transmitter that converts an excitation signal into an ultrasonic signal; an 
ultrasonic receiver positioned to receive the ultrasonic signal transmitted by the 
ultrasonic transmitter and that generates a received signal that is a nonlinear function 
of the excitation signal; an excitation signal generator in communication with the 
ultrasonic transmitter that generates an excitation signal that substantially maximizes 
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the signal generated by the ultrasonic receiver based on a specified constraint on the 
excitation signal; and a processing system in communication with the ultrasonic 
receiver for processing the signal generated by the ultrasonic receiver into an image of 
tissue disposed between the ultrasonic transmitter and the ultrasonic receiver. 

[0024] The signal generated by the excitation signal generator may be derived from a 
nonlinear model of the nonlinear function. The nonlinear model may be developed 
from a comparison of a transmitted test signal transmitted by the ultrasonic transmitter 
and a received test signal generated by the ultrasonic receiver. 

[0025] These processes and system may also be used in transmit-receive systems 
that operate outside of the ultrasonic range. 

[0026] These as well as still further features, benefits and objects will now become, 
clear upon an examination of the following detailed description of illustrative 
embodiments and the attached drawings. 

Brief Description of Drawings 

[0027] Fig. 1 is a block diagram of a transmit-receive system. 

[0028] Fig. 2 is a block diagram of a process for determining the optimum excitation 
signal for a transmitter in a non-linear transmit-receive system. 

[0029] Fig. 3 is a diagram of a Laguerre Volterra Network. 

[0030] Fig. 4 is a diagram of parallel cascades of linear filters and their associated 
nonlinearities. 

[0031] Fig. 5 is non-linear model of a transmit-receive system using a linear filter and 
a static nonlinearity. 

[0032] Figs. 6(a) - 6(d) are traces of measured and predicted signals in a non-linear 
transmit-receive system. 

[0033] Fig. 7 is a block diagram of a non-linear model of a non-linear system using 
principal dynamic modes. 

[0034] Figs. 8(a) - (d) are traces of optimal and pulse excitation signals and their 
corresponding received signals. 
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Detailed Description of Illustrative Embodiments 
[0035] Fig. 1 is a block diagram of a transmit-receive system. 

[0036] As shown in Fig. 1 , an excitation signal generator 101 may generate an 
excitation signal. This signal may be a steady-state signal, a series of pulses, or any 
other type of signal. 

[0037] The signal generated by the excitation signal generator 101 may be delivered 
to a transmitter 103. The transmitter 103 may be an ultrasonic transmitter or any other 
type of transmitter. The transmitter may be a single transmitter or an array of 
transmitters. In the case of an array of transmitters, separate signals may be delivered 
to each element in the array from the excitation signal generator 101. The separate 
signals may be the same or different. 

[0038] The signal transmitted by the transmitter 103 may travel through a 
propagation medium, such a breast suspended within acoustical coupling fluid, and 
then arrive at a receiver 105. The receiver 105 may be an ultrasonic receiver or any 
other type of receiver. As with the transmitter 103, the receiver 105 may be a single 
receiver or an array of receivers. 

[0039] The signal generated by the ultrasonic receiver 105 (or signals in the case of 
an array) may be directed to a processing system 107. The processing system may 
analyze the received signal, possibly along with the excitation signal, to extract 
information that is of interest. In one embodiment, the information of interest may be a 
three-dimensional image of a breast interposed between the transmitter 103 and the 
receiver 105. More details concerning transmission mode tomographic imaging for 
breast mammography are set forth in U.S. Patent Application Serial No. 10/11 7,763, 
filed April 5, 2002, entitled "High Resolution 3D Ultrasonic Transmission Imaging," U.S. 
Patent Application Serial No. 10/117,468, filed April 5, 2002, entitled "Nonlinear 
Processing for Mitigation of Diffraction Effects," U.S. Provisional Application Serial No. 
60/396,516, filed July 16, 2002, entitled "Support Bra for Ultrasonic Breast Scanner," 
U.S. Patent Application Serial No. 10/382,235, filed March 5, 2003, entitled "Multi- 
Mode Processing for Ultrasonic Imaging," U.S. Provisional Application Serial No. 
60/478,981, filed June 17, 2003, entitled "Targeted Ultrasonic Hyperthermia With A 
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High-Resolution 3-D Fixed Coordinate System," and U.S. Patent Application filed June 
27, 2003, entitled Three-Dimensional Ultrasonic Mammography Scanning," attorney 
docket number 64693-070. The entire content of all of these applications is 
incorporated herein by reference. 

[0040] Fig. 2 is a block diagram of a process for determining the optimum excitation 
signal for a transmitter in a nonlinear transmit-receive system. 

[0041] The transmitter 103 may be excited with a test signal, as reflected by an 
Excite Transmitter With Test Signal step 201. The test signal may be a wideband 
signal. It may be bandwidth-limited white noise, such as white noise with a bandwidth 
of 4 - 24 MHz. It may be a chirp. The wideband test signal may be selected to cover 
the bandwidth and dynamic range over which the transmitter 103, the receiver 105 and 
their associated electronics are configured to operate during normal use. The test 
signal may be any other form of signal. 

[0042] The test signal may travel through the anticipated propagation medium, such 
as acoustic coupling fluid containing a breast dangling therein. The test signal may 
then reach the receiver 105. The receiver 105 may then generate a received test 
signal, as reflected by a Receive Test Signal from Receiver step 203. 

[0043] The excitation signal and the signal generated by the receiver may be stored 
for later analysis. If desired, the signals may be digitized either before or after their 
storage. A time-bandwidth product of approximately 1 ,000 may be used for this 
purpose which, in the case of an excitation signal having a bandwidth of approximately 
4-24 Megahertz, may be approximately 100 MHz. 

[0044] The signal generated by the receiver 105 may be a nonlinear function of the 
excitation signal generated by the excitation signal generator 101. The nonlinearity 
may be caused by nonlinearity in the transmitter 103, the receiver 105, the excitation 
signal generator 101, the propagation medium, and/or by other areas of the system. 

[0045] A nonlinear model of the dynamic transformation between the excitation 
signal and the signal generated by the receiver may be developed, as reflected by a 
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Develop Nonlinear Model step 205. This may be a mathematical model that is derived 
using any type of nonlinear system identification method. 

[0046] Based on the nonlinear model that is developed, an optimum excitation signal 
may be determined that optimizes the signal generated by the receiver 105 based on a 
constraint on the excitation signal that is delivered to the transmitter 103. Any 
technique for determining the optimum excitation signal based on the nonlinear model 
may be used. • 

[0047] The constraint on the excitation signal may include a constraint on its peak 
amplitude, power and/or other characteristics. 

[0048] The excitation signal may be optimized to maximize the amplitude, power 
and/or other characteristics in the signal generated by the receiver 105. 

[0049] The optimum excitation signal may then be applied to the transmitter 103, as 
reflected by an Apply Optimum Excitation Signal step 209. To this end, the excitation 
signal generator 101 may be configured to generate the optimum excitation signal 
during the normal operation of the transmit-receive system. 

[0050] Fig. 3 is a diagram of a Laguerre-Volterra network (LVN). This is an example 
of a nonlinear system identification method that may be used as part of the process of 
developing a nonlinear model of the transmit-receive system. 

[0051] In Fig. 3, x(n) represents the excitation signal generated by the excitation 
signal generator 101 and y(n) represents a prediction of the signal that is generated by 
the receiver 105 as a result. 

[0052] The LVN method may combine expansions on the discrete-time Laguerre 
basis with feed-forward artificial neural networks (FANN) using polynomial activation 
functions. It may follow the standard architecture of a single-layer, fully connected 
FANN with three distinctive features: 

[0053] (1 ) The Laguerre-filter bank bj (/=0, 1 , . . . , L) that preprocesses the 

input, x(n); 

[0054] (2) The polynomial activation functions, ft, (/c=1 ,2, . . . , K) in the 

hidden units; and 
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[0055] 



(3) The non-weighted summative output unit. 



[0056] The architectural parameters of the LVN may consist of the number of 
Laguerre filters (LFs), hidden units (HUs) and the degree of the polynomial activation 
functions. The selection of these structural parameter values can be performed by 
successive trials in ascending order, moving from lower to higher numbers, using the 
normalized mean-squared error (NMSE) of the output prediction achieved by the 
model as a criterion. 

[0057] The determination of the optimal number of LFs for a specific system can be 
assisted by observing the relative magnitude of the Laguerre expansion coefficients. 
Likewise, the determination of the optimal number of HUs can be assisted by 
observing the relative magnitudes of the polynomial coefficients, when the inbound 
weights for each HU are normalized to a unity sum of squares. 

[0058] Once the LVN architecture is set, the training of the LVN may be performed 
with input-output experimental data to estimate the values of the unknown network 
parameters, including the weights w k j, the polynomial coefficients c qik , the Laguerre 
parameter a and the offset y 0 . The network training may involve minimizing the mean 
square error between the measured output and the LVN predicted output y(n) for the 
corresponding input. For effective training, the input-output data may be broadband 
and cover the entire bandwidth and dynamic range of the system. 

[0059] The output of the ;-th discrete-time Laguerre filter with impulse response 
function bj(m) may be given by the discrete convolution: 



[0060] 



Vj(n) = ^bj(m)x(n-m) 



(D 



[0061] where 



[0062] 



bj(m) = a 



(m-j)/2 



,=0 \l )\l J 



(2) 



[0063] and a is the discrete-time Laguerre parameter. 
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[0064] The input of the /c-th hidden unit may be a weighted sum of the Laguerre filter 
bank outputs: 

[0065] u k (n) = iX,v/«; (k=l,2,...K) (3) 

7=0 

[0066] The corresponding output of the k-ih hidden unit may be given by the 
polynomial activation function: 

[0067] h(n)=f d c vk utM (k = l,2,...K) (4) 

q=0 

[0068] The network output may be given by a non-weighted summation of the hidden 
unit outputs including a trainable offset value, y 0 : 

[0069] y(n) = j]Z k (n) + y 0 (5) 

[0070] The training of the network parameters may be performed iteratively using a 
gradient descent method for a quadratic cost function defined by the square of the 
output prediction error. 

[0071] As shown in Fig. 4, the LVN architecture shown in Fig. 3 may generally be 
expressed as K parallel cascades of linear filters Lu.Xk and their associated static 
nonlinearities Nf . . . N k . Fig. 4 shows the equivalent version of the LVN where the input 
x(n) goes through a set of K linear filters (determined by the weights w k j) and their 
associated nonlinearities (determined by the polynomial coefficients c qi k). The output 
may be formed by summing all of the outputs of the static nonlinearities. 

[0072] More details concerning the use of a Laguerre-Volterra network as a nonlinear 
system identification method are set forth in V.Z. Marmarelis, "Modeling Methodology 
for Nonlinear Physiological Systems," Annals of Biomedical Engineering, Vol. 25, pp. 
239-251 (1997), the entire content of which is incorporated herein by reference. 

[0073] Once the network parameters are found, the corresponding model of kernel 
functions may be constructed in accordance with well known techniques. The obtained 
mathematical model may describe the nonlinear dynamic transformation of the input 
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excitation signal into the output signal generated by the receiver in the transmit-receive 
system. 

[0074] In one experiment, the test signal was a wideband chirp covering a bandwidth 
from approximately 4 to 20 MHz over a dynamic range of ± 50 volts. The received 
signal was recorded over 5,000 samples at 20 nsec sampling intervals. 

[0075] The excitation waveform was generated from an arbitrary waveform 
generator, such as an HP 33250A from Hewlett-Packard in Palo Alto, California, with a 
pre-designed wideband waveform covering this bandwidth and dynamic range. The 
system output was measured at the receiving transducer and digitized using an 
oscilloscope, such as the TDS5054 Digital Phosphor Oscilloscope from Tektronix in 
Beaverton, Oregon. 

[0076] The digitized transmitted and received signals were analyzed using the LVN 
nonlinear system identification method described above to obtain a nonlinear model of 
the dynamic transformation of this transmit-receive system for any excitation waveform 
into the corresponding received waveform. 

[0077] Application of the LVN modeling methodology to the collected data yields the 
nonlinear model shown in Fig. 5. This is comprised of the excitation waveform 501 
being delivered to a linear filter 503, followed by a cubic static nonlinearity 505 and a 
received signal 507. 

[0078] With the LVN structural parameters of 1 7 Laguerre functions, 2 hidden units, 
and a third-degree polynomial activation function, the normalized mean squared-error 
of the model prediction with respect to the measured output data was 5.9% after 100 
iterations. In all trials with different structural parameters, only one hidden unit was 
found to be active. This allowed the model to be simplified to the single cascade 
shown in Fig. 5. 

[0079] Fig. 6 shows the measured receiver output and the LVN model predictions 
along with their power spectra. The accuracy of the LVN model prediction is evident in 
both the time and the frequency domains. 
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[0080] The optimal excitation waveform may be derived from this model on the basis 
of the matched-filter principle, that is, by determining the time-reversed impulse 
response of the linear filter of the cascade model, scaled according to an input power 
constraint. The optimal excitation pulse may match the frequency response of the 
linear filter. The static nonlinearity may determine the final gain of the output and 
indicate that doubling of the input more than doubles the output (supralinear relation). 

[0081] Fig. 7 is a block diagram of a nonlinear model of a nonlinear system using 
principal dynamic modes. The Pi .. . P k blocks represent the principal dynamic modes 
(PDM) of the system based on the kernel functions. The ui(t) . . . urft) represent a 
convolution of the input x(t) with each of the principal dynamic modes. The f block 
represents the static point-to-point mapping. 

[0082] Mathematically, y(t) in Fig. 7 may be expressed as: 

[0083] y(t) = f{u l (t),..,u k (tj\ (6) 

[0084] The input signal x(t) may be expressed as: 

[0085] *(0 = £vU-0 (?) 

[0086] where /\, are scalers. 

[0087] By applying the input signal, Uj(t) may be expressed as: 

[0088] Uj (f) = Z^PjWptT-WT (8) 

1=1 0 

[0089] =i>,.Z,-,,(r) (9) 

[0090] and y(t) may be expressed as: 

[0091] Xt) = f& 9 ~,*>i} (10) 

[0092] where Zjj (t) are known. 

[0093] y(t) may then be optimized over the parameters {X Jt ...,X k (subject to a 
constraint: 
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[0095] For example, the values (Ai A K ) may be determined for the maximum value 

of y(t) over t is greatest subject to the constraint: 



[0097] Other optimization criteria, such as maximizing the power of y(t) and/or 
constraints such as 0 < Ai < 1 , can also be applied. 

[0098] As indicated, these are examples of approaches for determining an optimum 
excitation signal based on an input constraint. 

[0099] Figs. 8(a) - (d) are traces of optimal and pulse excitation signals and their 
corresponding received signals. Fig. 8(a) is a trace of the optimal excitation signal that 
was determined based on the principals discussed above; Fig. 8(b) is a trace of the 
signal generated by the receiver of this optimal excitation signal; Fig. 8(c) is a trace of 
a non-optimal pulse that was commonly used to excite a transmitter; and Fig. 8(d) is a 
trace of the signal generated by the receiver of this non-optimal pulse signal. As can 
be seen by a comparison of the traces in Figs. 8(b) and (d), use of the optimal 
excitation signal increased the amplitude of the received signal by about 35 times, 
without any significant increase in the amplitude of the input signal. 

[00100] Although a test signal has been described as a mechanism for obtaining the 
nonlinear model of the system, other approaches could be used instead. In short, what 
has thus-far been described are merely examples of the components, steps, features 
and benefits of the invention. The invention is not limited to these, but solely to the 
subject matter delineated by the following claims and their equivalents. 



[0096] 



^(£,) 2 = constant 



(12) 
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